Baseline study for higher moments of net— charge multiphcity distribution at 

RHIC energies 
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Lattice QCD models predict the presence of a critical point in the QCD phase diagram where the 
first order phase transition between the hadron gas and Quark-Gluon Plasma ceases to exist. Higher 
moments of conserved quantities, such as net-charge, net-baryon number and net-strangeness, and 
their products, are proposed as sensitive probes for locating the critical point. The higher moments 
of net-charge multiplicity distributions have been studied for RHIC energies for Au-Au collisions at 
^sii7 = 7.7 to 200 GeV using three event generators, viz., HIJING, UrQMD and THERMINATOR- 
2. We discuss the effects of centrality selection, resonance production, as well as different particle 
species on the net-charge moments and their products. Finally, we present the model predictions 
on higher moments which set the baseline for measurements at RHIC. 

PACS numbers: 25.75.-q 



I. INTRODUCTION 

The Relativistic Hcavy-Ion Collider (RHIC) at 
Brookhaven National Laboratory is a dedicated fa- 
cility for the search and study of Quark-Gluon 
Plasma (QGP) and to explore the QCD (Quantum 
Chromodynamics) phase transition. The phase di- 
agram, plotted in terms of baryon chemical poten- 
tial (/is) and temperature (T) depicts the transition 
from QGP to a state of matter consisting of hadron 
gas (HG). Lattice QCD based models reveal that 
at vanishing fis, the transition from QGP to HG 
is a simple crossover whereas at large hb, the 
phase transition is of first order [l|-0|- Thus a criti- 
cal point in the QCD phase diagram is expected at 
the point where the first order transition ends. In re- 
cent years, probing the location of the critical point 
has been one of the primary goals of the physics pro- 
gram at RHIC. The QCD phase diagram is scanned 
by varying the center-of-mass energy of heavy-ion 
collisions at RHIC, thereby changing the fiB and T 
in a controlled manner. It has been predicted that 
the location of the critical point could be signaled 
by the non-monotonic behavior in several fluctua- 
tion observables as a function of the beam energy. 

Lattice and other QCD based models predict that 
higher moments of conserved quantities, such as net- 
charge, net~baryon number, net-strangeness, and 
net-isospin distributions, provide the most sensitive 
probes of the critical point. This is because of the 
fact that at the critical point, thermodynamic sus- 
ceptibilities and the correlation length (^) of the sys- 
tem diverge [1, 0, IS ■ The moments of the distribu- 
tions, mean (M), standard deviation (cr), skewness 
(S) and kurtosis (k), are related to the higher or- 
der thermodynamic susceptibilities of the conserved 
quantities and also to the correlation length. For 
example, the variance, skewness and kurtosis have 
been shown to be related to the ^^-^ and re- 



spectively In order to cancel the volume term 
in the susceptibilities, different combinations of the 
various moments arc constructed, some of these are: 
M /a"^ , Sa and Kcr^. These products, as a function of 
beam energy, are expected to show non-monotonic 
behaviour near the critical point. 

Experimental study of net-charge fluctuations 
may have some the advantages over other conserved 
quantities. Fluctuations of net-charges are related 
to the net-baryon and net-isospin fluctuations @. 
Due to the problem with detection of neutrons in 
heavy ion experiments, it has been shown that 
net-proton may not be a good proxy for the net- 
baryon [13, [ll|. On the other-hand, the net-charge 
takes into account all charged particles produced in 
the heavy ion collisions. Other analytical calcula- 
tion suggests that the net-charge higher moment 
analysis may be less affected by effect of acceptance 
as compared to that of net-proton. Recent lattice 
QCD model estimations [l3| add new motivation for 
higher moments analysis to extract the freeze-out 
parameter from the lattice QCD calculation. 

To map the QCD phase diagram and to locate the 
critical point, the beam energy scan program [l3 - [l6j 
has started in the year 2010 at RHIC. The program 
aims to scan the beam energy from -^snn = 7.7 to 
39 GeV including 62.4 and 200 GeV, corresponding 
to baryon chemical potentials from 450 MeV to 20 
MeV [13. 

In this study, the QCD based model (HI- 
JING [13 ), transport model (UrQMD [l^) and 
thermal model (THERMINATOR-2 ^) are used 
to study fluctuations of net-charge. To understand 
the various physics processes, like jet and mini-jet 
production by HIJING, hadronic re-scattering and 
various resonance decay by UrQMD and thermal 
production of particles by THERMINATOR-2 are 
used as background for this analysis. The above 
models are blind to the critical phenomena in the 
heavy ion collision physics. Thus, these models help 
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to understand the baseline for the higher moments 
of the net-charge distributions. 

Rest of the article is organized as follows. In the 
next section, we introduce the net-charge moments. 
In section III, we study the importance of central- 
ity selection and show how centrality bin width cor- 
rection can be used. In section IV, we discuss the 
contributions of different particle species on the net- 
charge moments. In section V, resonance effects are 
studied. Finally in section VI, we present the prod- 
ucts of moments from different event generators and 
compare to that of the hadron gas model. We con- 
clude with a summary and outlook for future studies 
in this direction. 



II. MOMENTS OF NET-CHARGE 
DISTRIBUTIONS 

In an experimental scenario, both the positive and 
negative charges are measured in a finite acceptance 
on an event-by-event basis. From these, the net- 
charge (difference between the positive and negative 
charges) distributions are obtained in a narrow win- 
dow in centrality. These distributions are to be an- 
alyzed to obtain various moments. The first four 
moments are expressed as. 
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where = — N- is the net charge and 6N = 
{N— < N >). The variance depicts the width of 
the distribution, the skewness is a measure of the 
asymmetry of the distribution, and kurtosis gives 
an indication of whether the distribution is picked 
at the center or not. 

Net-charge multiplicity distributions are sensitive 
to centrality of the collision. Centrality is normally 
expressed in terms of impact parameter of the col- 
lision or by percentage of total cross sections (0-5% 
being most central) or by the number of participat- 
ing nucleons (A^part)- Different centrality classes are 
defined within a narrow range of iVpart- The av- 
erage value of iVpart for different centrality class is 
estimated from Glauber model calculations. 

The beam energy dependence of the net-charge 
distributions for central collisions has been studied 
using different models by choosing the charged and 
identified particles selected within pseudo-rapidity 
range of |?7| — 0.5 and transverse momentum range, 
0.2 < PT < 2.0 GeV/c. In Fig. [Ha), this is shown 
for the UrQMD event generator, for four different 



center-of-mass energies. With the increase of beam 
energy, the total number of charged particle mul- 
tiplicity increases and particle to anti-particle ratio 
tends to unity. Thus the mean of the distribution 
get closer to zero and the distribution gets wider for 
increasing center-of-mass energy. 

Centrality dependence of the net-charge distribu- 
tions can be studied by plotting the distributions for 
various centralities. Fig. [IJb) shows the net-charge 
distributions for three different centrality classes for 
Au-Au collisions at y^sJiN — 39 GeV, obtained from 
the UrQMD model. It can be seen that the width of 
the net-charge distribution decreases as we go from 
central to peripheral events. Moreover, the mean 
of the distribution also shifts to higher values in 
going from central to peripheral events. The mo- 
ments of the distributions are sensitive to the width 
of the centrality window, which needs to be dealt 
with properly. This will be discussed in the follow- 
ing section. 



III. CENTRALITY BIN WIDTH EFFECT 
ON HIGHER MOMENTS 



Moments of the distributions are obtained from 
net-charge distributions for a given centrality class. 
A given centrality class is a collection of different 
events having different impact parameters. Colli- 
sion at the given impact parameter represents cor- 
responding volume of the collision geometry. Due 
to cvent-by-event estimations, distributions corre- 
sponding to several different impact parameters con- 
tribute with a given centrality window. This re- 
sults in additional fluctuations within each central- 
ity class. Due to this effect, the value of the mo- 
ments gets affected with increase in centrality bin 
width. This centrality bin width effect can be see 
in Fig. [5] where we plot the first four moments 
of net-charge distributions for Au-Au collisions at 
^SNN = 39 GeV. This is plotted for three centrality 
bin widths corresponding to 2.5%, 5%, and 10%. It 
shows first two moments, M and cr are not so much 
affected for different widths, but the S and k's val- 
ues are very much affected by these effect. The vari- 
ation for larger bin width as a function of centrality 
is quite large. Due to this effect, these moments in- 
crease with increase in centrality bin width. This ef- 
fect needs to be corrected for moments analysis. To 
correct this effect, the direct weighted method [2l| is 
used. In this method instead of direct estimation of 
higher moments in a given centrality class, following 
weighted averages are used: 
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FIG. 1; The net-charge muhiphcity distributions ob- 
tained from the UrQMD model for Au-Au colhsions as a 
function of (a) centrality, expressed in terms of percent- 
age or cross section for ^snn=39 GeV, and (b) coUision 
energy for y^s^=7.7, 19.6, 39 and 200 GeV. 



where the index i runs over each multiphcity bin, 
Xi represents various moments for the i-th bin, and 
Hi is the number of events in the z-th muhiphcity 
bin. Ui = N is the total number of events in the 
centrality bin. Fig|3] shows corrected moments by 
using this method of centrality bin width correction. 
No centrality bin width dependence is observed in 
each of three different centrality classes. For the 
rest of the article, all the plots are made by using 
the centrality bin width correction of the moments 
and the statistical error of the above higher moments 
arc estimated by the Delta theorem [23 |. 



IV. EFFECT OF PARTICLE SPECIES ON 
HIGHER MOMENTS 

Charged particles, produced in the collision, 
mostly comprise of 7r+, tt", A'+, K~ , p, and p. All 






100 200 300 

<N > 

^ part ' 



FIG. 2: Mean, standard deviation, skewness and kurto- 
sis for Au-Au collisions at -^snn=39 GeV as a function 
centrality, expressed in terms of number of participating 
nucleons. These values are obtained for the net-charge 
distributions with centrality bins of 2.5%, 5% and 10%. 
In this figure, centrality bin width effect is not taken care 
of. 



of these particle species contribute to the shape of 
the net-charge multiplicity distributions and hence 
to the moments of the distributions. To understand 
the contribution of different particles to the final net- 
charge higher moments, distributions for net-pion, 
nct-kaon and net-proton are made separately, tak- 
ing the particle identification at the event generator 
level. This has been studied for different centralities 
and different colliding energies. On the other hand, 
the leptonic contribution on the net-charge multi- 
plicity distribution has very negligible effect [1^ . 

In order to study the centrality dependence of 
the contribution from particle species for the mo- 
ments of net-charge distributions, we have used 
the UrQMD event generator for Au-Au collisions at 
•v/sNN=39 GeV. All types of particles are selected in 
the same pseudo-rapidity {\ri\ = 0.5) window and 
same transverse momentum ranges (0.2 < < 
2.0 GcV/c). Fig.glshows the M, a, S and k values of 
the net-charge, net-pion, net-proton and net-kaon 
distributions as a function of number of participant 
nucleus. We observe that except for net-pion, the 
mean of the other three distributions increase in go- 
ing from peripheral to central collisions. The net- 
proton distributions dominate the mean of the net- 
charge distributions. For nct-pions, the mean shifts 
towards negative value for central collisions. The 
width of the net-charge distribution is dominated by 
the net-pion distribution. Widths of the net-kaon 
and net-proton distributions are close together. The 
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FIG. 3; Mean, standard deviation, skewness and kurto- 
sis for Au-Au collisions at ^snn=39 GeV as a function 
centrality, after the bin width corrections are taken into 
account. The resulting moments of the distributions are 
close to each other for all three (2.5%, 5% and 10%) 
centrality bins. 
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FIG. 4: Mean, standard deviation, skewness and kurto- 
sis for Au-Au collisions at ^snn=39 GeV as a function 
centrality, expressed in terms of number of participating 
nucleons, for net-charge, net-pion, net-kaon and net- 
proton. 



skewness for net-proton is obviously very different 
compared to other three sets, because of the differ- 
ence between proton and anti-proton production in 
the collision. The skewness for net-charge is similar 
those of the net-pions. Similarly, the kurtosis values 
for net-charges arc also close to those of net-pions. 
Thus different moments of the net-charge distribu- 
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FIG. 5: Products of moments, Sa (upper panel) and K.a^ 
(lower panel), plotted with respect to average number 
of participating nucleons, for net-charge, net-pion, net- 
kaon, and net-proton distributions. 



tions are affected differently by net-pions, net-kaons 
and net-protons. 

In order to understand the contribution of differ- 
ent species on the product of moments, in Fig. [5] 
we show the products, Sa and Ka^ as a function of 
centrality for different particle species. The Sa val- 
ues do not show any centrality dependence for net- 
charge or any other particle species. The Sa values 
for net-charge are close to those of the net-pions, 
where are those of net-protons are much larger. The 
values of Ka^ for net-charges are dominated by those 
of the net -protons. A strong centrality dependence 
of the numbers is seen where Ka^ values for net- 
charge and net-protons decrease in going from pe- 
ripheral to central collisions. For net-proton as well 
as for net-kaon, the na^ values are very close to 
unity with no centrality dependence. 

It is known that the particle production mecha- 
nism for different particles depends on the center-of- 
mass energy of the collision. Thus as a function of 
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FIG. 6: Products of moments, Sa (upper panel) and na^ 
(lower panel), plotted with respect to Au-Au collisions at 
different colliding energies. The moments are obtained 
from net-charge, net-pion, net-kaon, net-proton distri- 
butions. 



FIG. 7: The product of the moments, Sa (upper panel) 
and K.a^ (lower panel), plotted with respect to average 
number of participating nucleons for two cases: with de- 
cay of all resonances and (ii) without resonance decays 
taken into account. 



the beam energy, the distributions of the net-charge, 
net-pion, net-kaon and net-proton will be different, 
thus affecting the higher moments and their prod- 
ucts. We have studied this in the UrQMD model, 
for top central (0-5%) collisions and function 
of beam energy, as shovifn in Fig. [6l For Sa, the 
energy dependence of net-proton is very prominent 
compared to net-kaon and net-charge, whereas no 
energy dependence is seen for net-pions. The Sa 
values for net-charge are close to net-pions for most 
cases. For kct^, no significant energy dependence is 
seen for net-proton or net-kaon. The energy depen- 
dence of net-charge is similar to that of the net- 
pions. Thus for both of the quantities shown in 
Fig. [51 energy dependence of net-charge is domi- 
nated by the net-pions. 



V. RESONANCE EFFECTS ON HIGHER 
MOMENTS 

In heavy-ion collisions, the abundance of reso- 
nance production may have its effect in net-charge 
higher moments in the finite pseudo-rapidity range. 
To study the effect of resonance decay on the net- 
charge higher moments analysis, we have used the 
thermal model like, THERMINATOR-2, for Au-Au 
collisions at ^snn =200 GeV. In this model, we are 
able to track the particles in order to know its origin. 
In the process, we can study the effects of the reso- 
nances such as, S, A++, p, </>, and w and their anti- 



particles. We have studied two cases: (i) with de- 
cay of the resonances, (ii) without resonance decays, 
where all the parent particles without any decay are 
considered. In Fig. [3 we present the results for Sa 
and KCT^ of net-charge distributions, as a function of 
number of participating nucleons, with and without 
resonance decays. It is observed that for both the 
quantities, the behavior of the data with and with- 
out resonance decay do not change. The values in 
fact increases with resonance decays turned on. 



VI. MODEL PREDICTIONS 

Distribution of the number of charged particles 
as well as the net-charges are sensitive to particle 
production mechanisms. Different models use differ- 
ent methods for particle production. The HIJING 
model treats the heavy-ion collisions as a superposi- 
tion of nucleon-nucleon collisions. It can be used to 
study the effect of jets and mini-jets on the produced 
particles. UrQMD, is a hadronic transport model 
including strings. It has been used successfully to 
describe stopping power and hadronic re-scattering. 
THERMINATOR-2 gives a good description for the 
thermal model of particle production. These mod- 
els do not contain any physics which is sensitive to 
the critical point. Thus these models can provide 
a baseline study for the physics analysis specific to 
QCD critical point. 

Net-charge distributions for the three models, HI- 
JING, UrQMD and THERMINATOR-2 are con- 
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FIG. 8: The (a) mean, (b) standard deviation, (c) skew- 
ness, and (d) kurtosis, of net-charge distributions, plot- 
ted with respect to average number of participating nu- 
cleons for Au-Au collisions at ^snn =200 GeV. The re- 
sults from various models, such as HIJING, UrQMD and 
THERMINATOR-2 are plotted. Expectations from cor- 
responding central limit theorem values are also super- 
imposed in terms of dashed lines. 



structed and analyzed to obtain the moments. Fig.|5] 
shows the M, u, S and K cLS 3t function of centrality 
for Au-Au collisions at ^snn=200 GeV for all the 
three models. The M and a values increase in going 
from peripheral to central collisions for all the three 
models, whereas S and k, values decrease with in- 
crease in collision centrality. The mean value of net- 
charge is lower in THERMINATOR-2 as compared 
to HIJING and UrQMD models, but the standard 
deviation values arc similar in all three cases. The 
centrality dependence of S for THERMINATOR-2 
is also weaker compared to other two models. 

The centrality evolution of the higher moments 
can be understood better by invoking the Central 
Limit Theorem (CLT) [HjHj, which gives the de- 
pendence of the moments on the number of partici- 
pating nucleons. These are: 
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Here A^part is the average number of participating 
nucleons. The centrality evolution of the higher mo- 
ments for all the three models follow the trend of the 
CLT in all the three models. 

Connection to lattice QCD calculations can be 
made by carefully selecting the products of the mo- 
ments, such as, M/cr^, Sa and kct^. In Fig. [9] we 
show the M/a^, Sa and kct^ for 0-5% centrality 
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FIG. 9: The collision energy dependence of the (a) 
M/a^, (b) Sa and na^ from net-charge distributions for 
top central (0-5%) collisions. The results are shown for 
HIJING, UrQMD, and THERMINATOR-2. The pre- 
dictions from the Hadron Resonance Gas model for the 
net-charge are plotted in all cases. 



for Au-Au collisions at y/sNN^7 .7 to 200 GeV for 
the three models. In addition, prediction from the 
Hadron Resonance Gas (HRG) [201 model are super- 
imposed on the plots. The M/a"^, and Sa values de- 
crease with increase in colliding energy. All the three 
models are close to each other in these cases. The 
values of na^ remain consistent for all energies ex- 
cept at ^/IWn = 200 GeV for HIJING and UrQMD. 
The HRG prediction shows larger values than the 
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three models in Mja^. and S'cr, but in ko^ . HRG 
model is very close to the three models except at 
200 GeV. In the HRG and THERMINATOR-2, the 
thermal equilibrium consideration is taken into ac- 
count for the particle production where as in the HI- 
JING and UrQMD particles productions are taken 
place by Lund string model. This could be the rea- 
son for the discrepancy between thermal (HRG and 
THERMINATOR-2), and other models like HIJING 
and UrQMD. 

VII. SUMMARY AND OUTLOOK 

The study of higher moments of net-charge mul- 
tiplicity distributions are proposed to be sensi- 
tive to the location of the critical point. In this 
manuscript, we have presented the higher moments 
of the net-charge distributions and their products 
for RHIC energies using HIJING, UrQMD and 
THERMINATOR-2 models. Net-charge distribu- 
tions inherently contain contributions from net- 
pion, net-baryon (net-proton) and net-strangeness 
(net-kaon). To understand the effect of individ- 
ual particles species on the net-charge distributions 
we have also studied net-pion, net-kaon and net- 
proton contributions. We observe that mean of 
the net-charge is dominated by net-protons whereas 



standard deviation, skewness and kurtosis are af- 
fected more by net-pion distributions. The reso- 
nance decay contributions to the net-charge distri- 
butions are studied by THERMINATOR-2 model. 
There is a small increase in ko^ and S'cr of the net- 
charge distribution values due to the resonance de- 
cay in our acceptance. The various moments by the 
three models arc compared for Au-Au collisions at 
y^NN=200 GeV and also compared with CLT ex- 
pectations. Finally, all three model predictions are 
compared with expectation from HRG model for var- 
ious colhding energies from ^smn=T-^ to 200 GeV. 
Model calculations are compared for MjcP' ^ na^ and 
Sa for high (0-5%) central collisions. These models 
show regular behavior of the products of moments. 
These studies set the baseline understanding for the 
current experimental program at RHIC for locating 
the QCD critical point. The effects of global charge 
conservation, finite acceptance, and system size de- 
pendence on the net-charge distributions will be fo- 
cused on a future study. 
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